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Abstract 

Following the wellknown spacetime decomposition technique as applied 
to (d+1) dimensions we write down the equations of magnetohydrodynamics 
(MHD) in a spatially flat generalised FRW universe. Assuming an equation 
of state for the background cosmic fluid we find solutions in turn for acous- 
tic waves and also for Alfven waves in a warm (cold) magnetised plasma. 
Interestingly the different plasma modes closely resemble the flat space coun- 
terparts except that here the field variables all redshift with their time due to 
the expansion of the background. It is observed that in the ultrarelativistic 
limit the field parameters all scale as the free photon. The situation changes 
in the prerelativistic limit where the frequencies change in a bizarre fashion 
depending on initial conditions. It is observed that for a fixed magnetic field 
in a particular medium the Alfven wave velocity decreases with the number of 
dimensions, being the maximum in the usual 4D. Further for a fixed dimension 
the velocity attenuation is more significant in dust compared to the radiation 
era. We also find that in an expanding background the Alfven wave propaga- 
tion is possible only in the high frequency range, determined by the strength 
of the external magnetic field, the mass density of the medium and also the 
dimensions of the spacetime. Further it is found that with expansion the 
cosmic magnetic field decays more sharply in higher dimensional cosmology, 
which is in line with observational demand. 
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1. Introduction 

The present investigation is a continuation of our earlier work [1] where we studied 
the propagation of an electromagnetic wave in a magnetised (un) plasma medium 
in an expanding universe. So we will be rather brief in mathematical details as 
also in other presentations. The work is primarily motivated by two factors. Al- 
though plasma is the most significant component of the cosmic fluid, barring some 
disjointed works the much sought after union between the plasma dynamics and 
cosmology has not been pursued so far to its logical end in a systematic way. Sec- 
ondly during the early phase of the universe before recombination the cosmic fluid 
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was presumably mostly in the form of a plasma whereas in the cosmological con- 
text the higher dimensional spacetime is particularly relevant in the early universe. 
When considering the cosmic evolution it is believed that during expansion at tem- 
perature (T > rUe) the matter existed mostly in the form of an electron positron 
plasma at ultra relativistic temperature during 10~^s < t < Is. This is followed by 
element formation so then we get essentially plasma of ionised hydrogen in thermal 
equilibrium with radiation. At this stage the energy density of radiation exceeded 
that of matter and the period loosely termed radiation dominated era. With further 
cooling (r < rUe) due to expansion recombination of ionised plasma took place so 
that matter got decoupled from radiation and we call it matter dominated era. On 
the other hand it can be shown that Einstein's equations generalised to higher di- 
mensions [2] do admit solutions where as time evolves the 'extra' dimensions shrink 
in size to be finally invisible with the current state of experimental techniques so 
the present visible world is effectively four dimensional. So we have thought it fit to 
study some of the plasma processes in an expanding background in the framework 
of higher dimensional spacetime. And this we have attempted to do in a system- 
atic way through the assumption of an equation of state (EOS) for the background 
fluid such that both the ultrarelativistic and prerelativistic era can be dealt with 
in a rather unified manner. We here tacitly assume that the plasma field is weak 
enough not to back react on the cosmic expansion, which is determined by the EOS 
obtaining at the particular instant. Moreover both the plasma field equations and 
the GR being highly nonlinear we are forced to consider a linearised plasma model 
for mathematical convenience. 

As commented earlier scant attention has been paid so far for the formulation 
of equations of magnetohydrodynamics (MHD) in curved spacetime, mathematical 
complexity being a strong deterrent in all the approaches. To our knowledge, Hol- 
comb and Tajima [3l H], Dettmann etal [5], Banerjee etal [U [7] are some of the 
authors making significant contributions in this regard. Like our previous work here 
also we adopt the (d-|-l) decomposition of spacetime as originally formulated by 
MacDonald and Throne [S] . In this approach there exists a set of preferred fiducial 
observers with respect to which the field quantities like mass density, velocity and 
electromagnetic field are all decomposed and measured. To illustrate the effects of 
curvature and expansion on MHD results we consider the spatially fiat FRW met- 
ric generalised to (d+1) dimensions as the simplest background available. At some 
stage we also assume the existence of an external magnetic field small enough not 
to disturb the isotropy of the FRW background i.e., we consider the region of space 
smaller than the horizon where the magnetic field may be taken to be coherent. In 
fact magnetic fields have been detected and measured in our galaxy and in our local 
group with the help of Zeeman splitting as also Faraday rotation measurements of 
linearly polarised radio waves from pulsar sources. We find that with expansion the 
cosmic magnetic field dissipates more efficiently in HD cosmology compared to the 
4D one. Given the fact that our cosmology is manifestly homogeneous and isotropic 
this rapid decay is desirable from the observational point of view. In a generalised 
FRW type of metric we systematically discussed the propagation of photons in 
empty space, in both cold and warm plasma with or without external magnetic field 
and find exact solutions of analogous equations of motion. Though not specifically 



stressed in the previous works by others we find that the different asymptotic modes 
closely resemble their flat space counterparts except by a conformal factor, which 
unlike the Newtonnian case makes the frequencies time dependent. Interestingly 
all the modes scale exactly as the free photon case except the matter dominated 
magnetised modes where the Alfven wave frequencies redshift in a bizarre fashion, 
depending on the mass density of the medium, strength of the magnetic field and the 
number of dimensions. We also find that the Alfven waves can either be evanescent 
or propagating depending on some initial conditions. It is found that the medium 
allows only high frequency waves to pass and with increase in dimensions the cutoff 
frequency also increases. 

We proceed to tackle our problem as follows: section 2 deals with the (d+1) 
decomposed GRMHD equations in the FRW setting while in section 3 we discuss the 
accoustic and Alfven waves and its characteristics. The work ends with a discussion 
in section 4. 

2. MHD equation in an expanding universe 

We extend here the (3+1) decomposition of GR as formulated by Arnowitt, Deser 
and Misner (ADM) [9J to a higher dimensional space time of (d+1) dimensions. As 
the split formalism has been extensively discussed and used in the literature [31 S] 
we skip the formalities and proceed to write down the final GRMHD equations for 
a spatially fiat generalized FRW metric 

ds^ = - [dx^ + dy"^ + dz^ + dil^l) (1) 

(n = 5, 6, 7, . . . , d ) 

where A = A{t) is the scale function. In a recent communication [1] we discussed 
in some detail the plasma phenomena in generalised FRW universe in the absence 
of an external magnetic field. We would like to complete our search by considering 
a magnetised plasma in this section to eventually find an expression for the Alfven 
wave. Let an ambient, large scale magnetic field [10] be present in the plasma 
medium and assuming a perfect fluid approximation we study the MHD equations 
in a curved background. 

We assume that the homogeneous magnetic fleld 'B' is generated prior to the time 
of recombination. Such a fleld could generated during electroweak phase transition 
[TTj.We write 

B = Bo{t) + Bi{x,t) (2) 

where Bi denotes the flrst order perturbation to the magnetic fleld. In an earlier 
work one of us [T2] explicitly found the solutions of the equations of motions for 
this spacetime via an equation of state, P = ap { P = pressure, p = energy density) 
such that the scale factor is 



(3) 



With the extrinsic curvature scalar defined as 
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we finally write down the curl equations of Maxwell [8] (see Mcdonald et al for 
3+1 ) split for more details ) generalized to ( d+1) dimensions as 



and 



A('^-i)(V xE) = -i^(A'^S) 



yl("-i)(V X = -^{A^E) + ^A" 
cot c 
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where the symbols have the usual significance. Again assuming a perfect conducting 
fluid for which the resistive force will be zero, we write 



E = -A- X B 

c 



(7) 



Now the equation of charge conservation is given by, 

-A'^V.j 



We assume that a perfect gas adiabatic EOS holds for a (d+1) dimensional case also 
such that P = ap^ where 7 is an adiabatic constant. Following Wienberg [13] one 



can easily show that here 7 comes out to be 7 = and the relativistic enthalpy, 



h = l + — + — 



(9) 



Further pressure gradient , VP = hC'l'W p where the velocity of sound is given by 
Cg = Now the linearised HD force balance equation is given by. 
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(10) 



For both radiation and matter dominated cases poh t and hCl ~ t ""^^"^ (po 
= mass density). From equation (6), neglecting displacement current, we get. 



t^'j = - (V X S) 



(11) 



Using the above relations we get from HD force balance equation, 
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Using equation (2) we get, x V x = But, 



and 



d 

r'^B X V X g^{t'"^B) = f^^^Bo X V X V X i>i X 5o (13) 



v|-(t"'^Pi) = ~r'^^hC^,poV{V.v,) (14) 



dt 

Using the above equation we get, 

od'^vi , , , , dvi 

t^^rir + (rid-f + 3nd - 2n - 4) t—^ 
dt^ at 

+ {2nd — n — 2){nd'-f + nd -\ — n — 3)vi 
+ t'^-^'^'^VAi X V X V X X 

- t'("'^^+"'^-"-^)C,',V(V.vi) = Min{d-i -d- l)e+^-^'"^vi X VAi (15) 

where, VAi = -^—^-^^ is a time invariant vector and the suffix 'z' denotes an 
arbitrary initial instant. 

2(l-da) 



It has not escaped our notice that the right hand side (RHS) of the equation (15) 
disappears whenever we consider an ultrarelativistic plasma in higher dimensional 
space time i.e., 7 = irrespective of the value of n, the expansion rate of the back 
ground FRW metric. The term Min{d'y — d — l)t"+^~^"''^fi x VAi remained unnoticed 
in the work of Holcomb et al. [3j but it was ffist noticed in the work of Sil et al [6] . 
This term will disappear for ultrarelativistic case only. 

As in Newtonian mechanics the relativistic Alfven wave velocity va is defined as 

VA = = VMi'-'^" = VMt-"^ (16) 

where VAi is the invariant Alfven wave velocity at any instant t = ti. For the 
prerecombination era in (d+l)-dimensional spacetime, a = ^ for perfect fluid and 

then Va ~ t~<*+i. Hence the Alfven velocity depends upon dimensions and decreases 
with it. It is maximum at 4D. So the Alfven velocity actually decreases in higher 
dimensions. 

On the other hand for a fixed the Alfven velocity varies as t^^ for the 
postrecombination era [a = 0). Alternatively, for the case of very large number of 
dimensions, the damping asymptotically reaches t~^, a form set for the postrecom- 
bination era also. It is to be noticed that the velocity for a = is independent 
of the number of dimensions, unlike the radiation case. The Alfven wave velocity 
decreases more sharply with time in dust (see figure 1). 



From equation (5) and (7) it follows that §-^{A^Bo) = ; or, ^(t^So) = 0. For 
radiation case {a = i), (-Bo)rad ~ with respect to comoving co-ordinates. 
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Figure 1: va ^ t graph 



It was pointed before that poh ~ both in higher dimensions and 4D and 
also in radiation and dust case. For the prerecombination era in (d+1)- dimensional 
spacetime Bo decreases with increasing dimensions, being maximum in 4D. Again 
for the postrecombination era (i.e, a = 0) the ambient magnetic field Bq ~ for 
a fixed 'd'. Bq reaches asymptotically for a very large number of dimensions in 
radiation case and the nature is same as the postrecombination era (see figure 2). 




Figure 2: Bq ^ t graph 



The fact that va decreases in an expanding background owes an explanation as 
follows : we know from classical idea that Alfven wave is generated as in a transverse 
vibration in a string as an interplay between the unbalanced tension and inertia. 
When looked upon as lines of force the tension in the string comes out to be 
the energy density of the magnetic field. Since with higher dimensions the volume 
increases the energy density, and consequently the magnetic field B also gets diluted. 
So the Alfven velocity also decreases with increasing dimensions. 



One can look at from a different standpoint also. As a consequence of flux 
conservation we may write 



dt J dt J 



e {e.m.f.) 



(17) 



Since the uniform ambient magnetic field does not produce an e.m.f we see that 

2 

if e = 0, as the magnetic field must decay in time t i+" with respect to co-moving 
co-ordinates. This result is also same for 4D well as for higher dimensions. 

Thus as a consequence of magnetic flux conservation the field lines in an expanding 
universe are conformally diluted. 

The theoretical result that with expansion the cosmic magnetic field decays much 
faster in HD compared to 4D may have non trivial implications in astrophysics. In a 
recent communication Gang Chen etal. p3j argued that if a primordial magnetic field 
did exist it should create vorticity and alfven waves and finally cosmic microwave 
background (CMB) anisotropics. They later went to discuss the constrints that can 
be placed on the strength of such a field with the help of the CMB anisotropy data 
from WMAP experiment. The origin of the observed magnetic galactic field of a 
few /iG apparently coherent over a 10 kpc scale continues to evade satisfactory the- 
oretical explanation. A plausible explanation may be the consequence of nonlinear 
amplification of a tiny seed field by galactic dynamo process [15]. CMB anisotropy 
data apparently suggests a seed field of nG [H] strength. Now the dynamo mech- 
anism being so efficient vis-a-vis amplification of the seed field at the early era it is 
quite likely to shoot up the seed field beyond the current value of a few /iC However 
that may introduce considerable anisotropy [16] in the manifestly homogeneous and 
isotropic FRW model. And here the HD model makes its presence felt in allowing 
the amplified field to decay much faster. 

Now, putting 7 = and n = -^^^ in equation (15) we get finally, 



The above equation is quite complex. It can be simplified, however, if we consider 
the interesting case of the wave moving perpendicularly to both magnetic field and 
wave vector. In this case 
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and 
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so the above equation gives 



id - Ma - 2 dvi 
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(21) 



where the velocities, C^j originates from the acoustic pressure and vai from the mag- 
netic field. This is a Bessel type of equation of order, p = {a ^ 0). 
The equation (21) yields a general solution 



-3d+5da + 2 
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(22) 



Special Cases: We now briefly discuss the equation (22) for the special cases of 
ultra relativistic (T >> rric) and also stiff fluid back ground. Evidently the prerela- 
tivistic case, a = is not obtainable from the above equation and needs a separate 
treatment. So we take radiation and for the sake of completeness the stiff fluid case 
also. 



I. Radiation Case : Putting a = ^ the equation (22) reduces to 



(3d- 7) 



Vl = t W+^H, 



ipi) 



d -\- 1 , r, , , 2 



(23) 



where the order of the Henkel coefficient is given by pi = the frequency of 

the oscillation ~ t~d+T and {v\^ + C^)^ is known as phase velocity. This mode 
corresponds to the magnetosonic or fast Alfven mode because both velocities are 
simultaneously present. The result is identical to the fiat space case also except 
that frequency redshifts here and is model dependent. 



This wave is somewhat similar to an electromagnetic wave, since the time varying 
magnetic field is perpendicular to the direction of propagation but parallel to the 
magnetostatic field, whereas the time varying electric field is perpendicular to both 
the direction of propagation and the magnetostatic field which is a longitudinal 
wave. However, since the velocity of mass flow and also fluctuating mass density 
associated with the wave motion are both in the wave propagation direction this 
wave is called the magnetosonic longitudinal wave. The phase velocity {v\^ + 0^)^ 
is independent of frequency so that it is a non dispersive wave. We see that the fluid 
velocity depends upon dimension also. It will reduce to Holcomb et al [3] work for 
d = 3. As noted earlier the attenuation of the wave propagation is more prominent 

/ -3d+5 \ 

in the higher dimensional regime compared to the 4D easel fi ~ t2{d+i) j \^ ^j^^ 



dimension 'd' increases indefinitely the damping finally becomes independent of 'd', 

3 

d — >■ cxo and vi t~2 (see figure 3). Now if one can consider a relatively cold plasma 
such that the 'pressure' due to magnetic field far exceeds the acoustic pressure one 
may take Csi = 0, this mode is known as pure or shear Alfven wave. In this case only 
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Figure 3: ui ~ t graph 



magnetic field is present, but the nature of wave will be same as magnetosonic wave. 
Another interesting feature comes out of the analysis. We have shown earlier pLj like 
others in the related field that a free photon redshifts in the ultrarelativistic case and 
the field variables in the different plasma models also redshift identically. We get 
exactly similar result for the Alfven mode also. Though not specifically emphasized 
in the earlier works this result, in our opinion, is not model independent. The 
fact that a conformally flat FRW mechanism is chosen in al our works all the field 
variables here scale with the conformal scale factor. A less symmetrical spacetime 
would possibly yield different results. We shall, however, subsequently see that the 
situation drastically changes in the non relativistic background. 

On the other hand for very low magnetic field the acoustic pressure dominates 
and we get pure acoustic wave. 

II. Stiff Fluid [a = 1): For the sake of completeness we also consider the stiff 
fluid case characterised by a = 1. Equation (22) gives {Csi= c = velocity of light ). 

vi=t'^H^^^ [{v\i + c^)kit] (24) 

However, the stiff fluid being not much physically relevant we skip discussion of 
this result. 



III. Dust Case (a = 0): Here equation (21) becomes, 



orfV 4c/ -2 dvi [2(^-1) . , 1 

* V + —'IT + + = ° <2^) 

It is clear that the form of the two equations corresponding to radiation and dust 
cases are distinctly different - one resembling the Bessel form while the other to an 
Euler type. In this case if the magnetic field is very strong so that the magnetic 
pressure is very much larger than the fluid pressure, then phase velocity of the 
magnetosonic wave becomes equal to the Alfven wave velocity vai- The second one 
has three solutions depending upon the discriminant D = {d — 2)^ — {2dkfv\^: 



(a) For D>0, 

V, = Cir^[(3o;-2)+7D] ^ ^^^-^[(3d-2)-VD] ^26) 

where ci and C2 are constants. This is a general solution of equation (21). Here the 
above mode always dies out without much propagation. 



(b) For D = Q, 



(3d-2) 

Vi = {C3 + C4\nt)t 53- (27) 



C3 and C4 are constants. It is apparent that fiis least in the usual 3D. This result is 
somewhat intriguing when posited against the fact that the phenomena of damping 
is due to the background expansion, which is maximum in 3D. So other process may 
be in work besides the expansion rate in this case. 

(c) For D <0, 

3d-2 I y U ^ \ 3d-2 , V ^ , 1 / n 

vi = c^t 2d COS -T^-lnt + Cat 2d sin -—\nt (28) 





where C5 and Cg are constants. So we get damped oscillatory mode here, although 
the periodicity is not in cosmic time but in its logarithm. Hence the signature of 
of D which, in turn, depends on the dimension, external magnetic field and the 
density of the medium is crucial in determining whether a particular mode will be 
evanescent or propagating. In fact for oscillation to take place one should have 
kiVAi > Now for dust VAi is independent of 'd'. So as number of dimension 

increase only very high frequency waves can propagate through a particular plasma 
medium. The extra dimensions, so to say, try to filter the frequency range. This 
finding is interesting but its astrophysical implications are yet to be looked into. The 
shear Alfven waves are stationary in a certain frequency range which depends upon 
D. For higher frequencies, the waves do propagate; however, they are not oscillatory 
in the cosmic time, but in its logarithm. 



3. Discussion 



As physics in higher dimension is increasingly becoming an important and at the 
same time a distinct area of activity we are primarily interested in this work to 
investigate the effect of higher dimensions in the propagation of an electromag- 
netic wave in plasma media with an expanding background. As commented earlier 
in the introduction this work may be looked upon as continuation of our earlier 
work where acoustic and Alfven modes are not considered. Here the background 
cosmology chosen is the simplest possible - FRW spacetime generalized to (d+l) 
dimensions. So the results we obtain difi^er essentially in some quantitative aspects 
without introducing much qualitative differences. 



The interesting results may be briefly summarised as foUows: We know that the 
expanding background generates a sort of damping such that the plasma variables 
like the amplitude of the wave, strength of the magnetic field etc. decay. Unlike 
mechanical damping the expansion creates a thinning of density of magnetic lines 
of force so that the velocity of the different modes keep on decreasing with time. 
Naturally for a contracting mode just the opposite will happen. We observe that 
this damping is the least in the usual 3D space compared to the HD space. As 
commented in the introduction this exercise may be looked upon as continuation 
of our earlier work where acoustic or Alfvcn waves are not discussed. The locally 
measured frequencies resemble that of flat spacetime, though due to the expansion 
of the universe they do decay in time. We have found like others in the related held 
that in the ultra relativistic limit all modes redshift in the same scale as free photon. 
In the non relativistic limit, however, the situation drastically changes. Different 
modes redshift at different rates depending on the plasma density, magnitude of 
the external magnetic fleld and interestingly on number of dimensions. So here 
dimensions may embed a deflnite signature. These factors also dictate whether a 
particular mode becomes evanescent or can propagate through the medium. We flnd 
that as dimension increases only very high frequency waves can propagate but they 
are not oscillatory in cosmic time but in its logarithm. Finally coming back to the 
vexed question of a primordial magnetic fleld we flnd the fleld decays much faster 
in HD spacetime, which is a good news for the observed isotropy and homogeneity 
of the present universe because large magnetic fleld introduces anisotropy. 

To end a flnal comment may be in order. We know that at the early phase when 
consideration of both higher dimensions and plasma phenomena are interesting the 
fluid was far from being perfect. Viscosity and heat processes were signiflcant. Pos- 
sible extension to this work could include formulating these equations for imperfect 
MHD including such effects as viscosity and flnite conductivity. Nonlinear modes 
such as shocks could also be investigated. Judging by the complexity of the solutions 
it seems likely that more complicated modes may not be analytically solvable and 
may require extensive numerical analysis. Our simple analysis has kept intact the 
basic essence of MHD theory, highlighting the similarities and occasional differences 
with standard flatspace MHD, which complicated numerical analysis may mask. 
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